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Abstract 

In  the  previous  work  [Metzler  et  al.,  Phys.  Plasmas  6,  3283  (1999)]  it  was  shown 
that  a  tailored  density  profile  could  be  very  effective  in  smoothing  out  the  laser  beam 
non-uniformities  imprinted  into  a  laser-accelerated  target.  However,  a  target  with  a 
smoothly  graded  density  is  difficult  to  manufacture.  A  method  of  dynamically  producing 
a  graded  density  profile  with  a  short  shaping  laser  pulse  irradiating  a  foam  layer  on  top 
of  the  payload  prior  to  the  drive  pulse  is  proposed.  It  is  demonstrated  that  the  intensity 
and  the  duration  of  the  shaping  pulse,  the  time  interval  between  the  shaping  pulse  and  the 
drive  pulse,  and  the  density  ratio  between  the  foam  and  the  payload  can  be  selected  so 
that  the  laser  imprint  of  the  drive  pulse  is  considerably  suppressed  without  increasing  the 
entropy  of  the  payload.  The  use  of  the  foam-plastic  target  and  a  shaping  pulse  reduces  the 
imprinted  mass  perturbation  amplitude  by  more  than  an  order  of  magnitude  compared  to 
a  solid  plastic  target.  The  requirements  to  the  smoothing  of  the  drive  and  shaping  laser 
beams  and  to  the  surface  finish  of  the  foam-plastic  sandwich  target  are  discussed. 
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I.  INTRODUCTION. 


Reduction  of  the  laser  imprint  remains  one  of  the  important  problems  of  the 
direct-drive  laser  fusion.  Even  with  all  the  improvements  of  the  optical  smoothing 
techniques,  the  residual  laser  beam  non-uniformities  imprinted  into  the  imploding  target 
could  be  amplified  by  the  Rayleigh-Taylor  (RT)  instability  to  amplitudes  large  enough  to 
prevent  high  gain  or  even  ignition. 12 

Most  approaches  to  making  advanced  imprint-resistant  targets  investigated  so  far 
are  based  on  using  either  low-density  foam  layers;3’4  external  x  rays  for  creating  a 
plasma  before  the  drive  laser  pulse,5’ 6  or  x  rays  from  high-Z  materials  in  the  target  itself 
heated  by  the  laser  beam  used  to  drive  the  target,2  or  some  hybrid  combinations  (such  as 
foams  preheated  by  x  rays4' 610  or  a  high-Z  layer  on  the  plastic  membrane  of  cryogenic 
deuterium-loaded-foam  targets11).  The  hybrid  approach  4  610  uses  x  rays  early  in  the 
pulse  to  pre-heat  the  outer  layers  of  the  target,  while  most  of  the  target  mass  is 
maintained  on  the  low  adiabat.  Alternatively,  a  short  laser  pre-pulse  before  the  arrival  of 
the  drive  pulse  could  be  used  for  pre-heating.12"15  The  purpose  of  pre-heating  in  all  these 
cases  is  to  decrease  the  density  and  to  increase  the  temperature  of  the  outer  layers  of  the 
target  some  time  before  the  drive  laser  pulse  arrives.  Both  factors  reduce  the  RT  growth 
rate:  the  lower  density  contributes  by  increasing  the  mass  ablation  velocity,  whereas  the 
higher  temperature  increases  the  smoothing  due  to  thermal  conductivity.  The  pre-heating 
must  be  confined  to  the  outer  layers  of  the  target  so  that  the  imploded  fuel  remains  on  a 
low  adiabat  consistent  with  ignition  or  high  gain.  One  should  also  mention  a  related  early 
idea  of  replacing  a  single  shaped  laser  pulse  with  a  train  of  shorter  pulses  to  improve  the 
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target  performance16  (the  concept  that  has  recently  re-emerged  as  the  picket  fence 
approach  to  shaping  the  low-intensity  foot  of  the  laser  pulse17). 

A  different  approach  to  the  mitigation  of  laser  imprint  is  the  use  of  tailored 
density  profiles,  first  suggested  in  Ref.  3.  The  method  of  density  tailoring  was  originally 
proposed  to  mitigate  the  RT  instability  in  Z -pinch  implosions.18  The  stabilizing  effect  is 
obtained  by  shock  waves  slowing  down  as  they  propagate  into  increasing  density.  The 
effective  gravity  near  the  shock  front  then  has  the  same  direction  as  the  density  gradient, 
which  makes  the  mass  perturbations  in  the  shocked  plasma  oscillate  rather  than  grow. 
Density  tailoring  seems  to  improve  radiative  performance  of  Z-pinch  plasma  radiation 
sources:  For  example,  the  cross-section  of  the  gas  jet  formed  by  a  double-shell  nozzle19 
(which  corresponded  to  an  appropriately  tailored  radial  density  profile)  was  shown  to 

give  the  most  significant  contribution  to  the  record-high  argon  K-shell  yield  (270  kJ) 

20 

recently  obtained  on  the  Z  facility  in  Sandia. 

Even  if  the  laser  target  initially  has  no  density  gradient,  the  mass  variation 
amplitude  will  oscillate  during  the  shock  transit,  as  predicted  by  theory21  and 
simulations,22  and  recently  confirmed  by  experiments  on  the  Nike  laser  at  the  Naval 
Research  Laboratory  (NRL). 23-24  The  additional  effective  gravity  provided  by  the 
density  tailoring  helps  increase  the  frequency  and  reduce  the  amplitude  of  such 
oscillations,  which  results  in  a  very  efficient  suppression  of  the  laser  imprint.  Simulations 
of  Ref.  3  demonstrated  that  a  planar  target  with  a  tailored  density  profile  is  much  more 
imprint-resistant  than  either  a  low-density  target  or  a  foam-plastic  sandwich  target  of  the 
same  total  mass  and  thickness.  There  is  some  experimental  evidence25  that  a  target 
approximating  a  graded  profile  with  a  three-stepped  density  structure  (layers  of  foam 
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with  densities  0.04,  0.12,  0.37  g/cm3,  each  5  jim  thick)  was  more  imprint-resistant  in  a 
wide  range  of  wavelengths  than  a  single  slab  of  a  low-density  foam. 

At  the  moment,  no  technology  is  available  to  manufacture  a  foam  target  whose 
density  is  smoothly  graded  in  one  (radial)  direction.  However,  the  same  physical 
mechanism  of  laser  imprint  mitigation  through  inverted  acceleration  could  be  used  if 
the  tailored  density  profile  were  produced  on  the  fly  dynamically;  by  a  separate  short 
laser  pulse  incident  upon  the  appropriately  designed  target  prior  to  the  drive  laser  pulse. 
Such  a  target  might  consist  of  a  sandwich  of  a  low-density  foam  on  top  of  a  higher- 
density  payload.  We  will  call  this  short  pulse  a  shaping  pulse  to  emphasize  its  purpose 
(shaping  the  density  profile),  differentiating  it  from  a  pre-pulse12"14  or  a  high-intensity 
spike,  a  picket  at  the  start  of  the  foot  preceding  the  main  drive  pulse.15  The  shock 
wave(s)  launched  by  the  drive  laser  pulse  then  propagate(s)  through  the  graded  density 
profile  created  by  the  expansion  wave  which  follows  the  first  shock  wave  produced  by 
the  shaping  pulse.  The  intensity  and  the  duration  of  the  shaping  pulse,  the  time  interval 
between  the  shaping  pulse  and  the  drive  pulse,  and  the  density  ratio  between  the  foam 
and  the  payload  can  be  selected  self-consistently  to  provide  a  shock  impedance  matching 
between  the  foam  and  the  payload.  The  shock  waves  generated  by  the  shaping  pulse  and 
the  drive  pulse  are  timed  to  merge  when  they  reach  the  contact  interface  between  the 
shock-compressed  foam  and  the  payload,  and  there  is  no  shock  reflection  from  this 
interface  because  the  density  of  the  shock-compressed  foam  is  close  to  the  density  of  the 
payload.  It  should  be  noted  that  the  purpose  of  the  shaping  laser  pulse  needed  for  this 
method  is  different  from  that  discussed  in  previous  works. 1217  Our  goal  is  not  primarily 
to  pre-heat  and  soften  the  outer  layers  of  the  target.  Rather,  we  intend  to  dynamically 
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generate  a  tailored  density  profile  that  satisfies  the  conditions  of  shock  impedance 
matching.  This  requires  a  more  complex  design  —  in  particular,  the  outer  foam  layer 
must  be  appropriately  designed  to  produce  the  tailored  density  profile.  This  design, 
however,  leads  to  a  much  higher  efficiency  of  imprint  mitigation  compared  to  that 
achieved  by  pre-heating  the  outer  layers  of  a  solid  target. 12-15 

This  paper  is  structured  as  follows.  In  Section  II  we  describe  how  the  shaping 
pulse  and  the  target  parameters  are  chosen  to  produce  the  desired  tailored  density  profile, 
achieving  at  the  same  time  the  shock  impedance  matching  at  the  foam/payload  interface. 
In  Section  III  we  demonstrate  that  the  target  designed  this  way  is  imprint-resistant  and 
find  the  limitations  imposed  on  the  non-uniformity  of  the  smoothed  shaping  pulse  and 
drive  laser  pulse  beams,  and  on  the  target  surface  finish.  In  Section  IV  we  conclude  with 
a  discussion. 

II.  TAILORING  A  DENSITY  PROFILE  WITH  A  SHAPING  LASER  PULSE 

This  Section  investigates  the  tailoring  of  a  density  profile  with  a  shaping  laser 
pulse  in  a  planar  foam-plastic  layer  (  sandwich  )  target  for  laser  imprint  mitigation.  The 
goal  is  to  find  the  conditions  under  which  this  concept  could  be  experimentally  tested  in 
planar  geometry,  similar  to  earlier  experiments.”’  ’  “  The  planar  target  emulates  the 
early-time  performance  of  a  small  part  of  a  spherical  target.  This  is  adequate  for 
modeling  the  imprint  mitigation,  since  the  perturbations  imprinted  into  the  accelerated 
plasma  from  the  optically  smoothed  laser  beam  are  mostly  in  the  short  wavelength  range 
and  the  target  moves  only  a  small  fraction  of  the  radius.  Nike  experiments  (laser 
wavelength  A  =  0.248  pm,  focal  length/=  600  cm,  diameter  of  the  laser  aperture  DA  =  15 
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cm),  are  typically  conducted  with  ISI  smoothing  at  1  THz  bandwidth,  which  puts  most  of 
the  laser  imprint  in  the  wavelength  range  20-30  pm  (see  Ref.  26),  whereas  the  flat  central 
region  of  the  Nike  beam27  is  400  pm  in  diameter  compared  to  a  typical  ICF  pellet  radius 
of  over  a  millimeter. 

What  is  referred  to  as  the  drive  laser  pulse  is  the  Nike  pulse  with  or  without  a 
low-intensity  foot,  which  sets  the  accelerated  payload  at  the  proper  adiabat  (for  details, 
see  Ref.  28).  The  energy  of  the  shaping  pulse  is  very  small  compared  to  the  energy  of  the 
drive  pulse. 

The  graded  density  profile  is  formed  in  a  process  called  impulsive  loading.  This 
can  be  idealized  as  a  fast  release  of  a  finite  energy  in  a  thin  layer  near  the  surface  x  =  0  of 

no  . 

a  half-space  x  <  0  filled  with  a  uniform  material  at  rest.  Impulsive  loading  has  been 
studied  since  the  1950s  when  it  was  visualized  as  an  explosion  of  a  thin  layer  of 
explosive  on  a  surface,  or  an  impact  of  a  thin,  light  plate  carrying  a  finite  kinetic  energy. 

A  short-duration  shaping  laser  pulse  delivering  a  finite  energy  to  a  thin  surface  layer  of 
the  foam  also  acts  as  an  impulsive  load.  The  impulsive  loading  produces  a  shock  wave 
that  propagates  into  the  target  in  the  negative  x  direction,  and  is  immediately  followed  by 
an  expansion  wave  which  gradually  reduces  the  shock  strength,  accelerating  the  fluid  in 
the  positive  x  direction.  If  the  shock  wave  is  sufficiently  strong,  and  the  ideal  gas  model 


is  a  good  approximation  of  the  equation  of  state,  then  the  flow  of  the  shocked  fluid  is 
self-similar.29,30  The  position  of  the  shock  front  at  time  t  is  given  by 
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where  /  and  x  are  the  intensity  and  duration  of  the  shaping  pulse,  respectively, 

Ep  =  I  x  is  the  energy  per  unit  area  delivered  in  the  impulsive  loading,  a  is  the  self¬ 
similarity  exponent,  <§„  is  the  dimensionless  factor  depending  on  the  equation  of  state. 
The  scaling  of  xs  (t)  with  E p  and  p0  is  the  same  as  in  the  case  of  a  planar  blast  wave/1 
In  the  latter  case,  the  dependence  of  xs  (7)  on  all  the  relevant  parameters  is  obtained 


directly  from  the  dimensional  considerations: 


(2) 


where  the  constant  corresponding  to  the  blast  wave  solution  of  Ref.  31,  is,  of  course, 
different  from  <§0  in  (1).  Proportionality  between  xs  (/)  and  (Ep  /  p0)'  in  (1),  of  course, 


follows  from  the  same  dimensional  considerations.  The  time  dependence  in  (1),  however, 
cannot  be  found  this  way,  since  the  self-similar  solution  under  consideration  belongs  to 
the  so-called  second  kind.29  The  self-similarity  exponent  a  is,  in  a  general  case,  an 
irrational  number  determined  from  the  condition  of  analyticity  of  the  self-similar  solution 
on  the  limiting  characteristic.  Note  that  the  shock  motion  in  the  limit  x  — >  0,  I p  — >  °o  , 

finite  E  ,  is  not  determined  by  the  energy  E p  instantly  released,  as  is  the  case  for  the 

blast  wave  (2),  but  rather  by  some  combination  of  E  and  x.  Obviously,  a  must  be  less 

than  one,  so  that  the  shock  wave  slows  down  as  it  propagates.  A  more  detailed  analysis  in 

Ref.  30  yields  the  following  inequalities  for  a: 

1  2  ™ 

—  <  cc  <  — .  (3) 

2  3 
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Comparing  (1)  and  (2),  we  see  why  a  should  not  exceed  2/3:  this  value  of  the  self¬ 


similarity  exponent  corresponds  to  the  case  of  a  blast  wave,  where  the  shocked  fluid  does 
not  expand  beyond  the  plane  of  symmetry,  x  =  0,  whereas  our  case  of  impulsive  loading 
permits  unlimited  expansion  in  the  positive  x  direction.29  Therefore  the  blast  wave  is 
supported  by  a  higher  pressure,  and  accordingly,  its  asymptotic  decay  is  slower. 

The  self-similar  density  profile  is  expressed  as 


P(x,  t)  — 


p0>  x<xs( 0; 

^|p0/V(r7),  x>xs(0, 

7-1 


(4) 


where 


(5) 


is  the  self-similar  coordinate,  and  the  function  N(rj)  satisfying  7V(T)  =  1  determines  the 
shape  of  the  profile.  The  time  dependence  of  the  self-similar  profile  is  thus  reduced  by  its 
stretching  along  the  horizontal  axis.  The  density  peaks  near  the  shock  front  and  decays 
with  increased  distance  from  it,  asymptotically  as  N(rj )  ~|  q  |“1/(1_a)  at  77  — >  .  There  is 

an  exceptional  case  of  7=  7/5,  for  which  the  self-similar  solution  is  found  explicitly:30 
a  =  3/5,  N(ri)  =  (5-4q)5'2,  -oo<q<\.  (6) 

The  dynamically  formed  density  profile  (4  )  satisfies  the  requirements  of  Ref.  3  to 
provide  an  efficient  imprint  stabilization.  Of  course,  the  impulsive  loading  is  not  fully 
equivalent  to  manufacturing  a  target  with  a  graded  density,  since  the  fluid  pressure  and 
velocity  are  nonzero  now,  also  featuring  the  self-similar  profiles.  However,  the  shock 
wave  launched  by  drive  laser  pulse  will  be  much  stronger  than  one  generated  by  the 
shaping  pulse.  Therefore,  the  pressure  and  velocity  variation  in  front  of  this  strong  shock 


wave  are  of  minor  importance,  and  its  propagation  is  governed  mostly  by  the  pressure 
driving  it  and  the  density  profile  ahead  of  its  front,  just  as  in  Ref.  3.  In  particular,  the 
shock  generated  by  the  drive  pulse  sets  the  payload  at  almost  the  same  adiabat  as  it  would 
without  the  shaping  pulse. 

To  provide  the  approximate  shock  impedance  matching,  the  peak  density  in  the 
graded  profile,  (y  +  l)p0  /(y  - 1) ,  must  be  equal  to  the  density  of  the  payload,  p, ,  i.  e. 


7~  1 

Po  =iPi- 

y+1 


(7) 


The  density  tailoring  can  only  help  stabilize  the  wavelengths  that  are  not  much  longer 
than  the  thickness  of  the  graded  density  layer.  The  initial  foam  layer  is  shock-compressed 
approximately  by  a  factor  of  (y  + 1)  /(y  - 1)  ,  hence  its  initial  thickness  L0  is  chosen  from 
the  condition 


L 


o 


(8) 


where  A  =  Am  is  the  longest  wavelength  to  be  stabilized.  For  instance,  assuming  that  the 
payload  is  made  of  a  solid  plastic,  p{  -  1.07  g/cm,3  Xm  =  30  pin  as  above,  and  7  =  5/3 
for  a  plastic  foam,  we  find  the  approximate  parameters  of  the  foam  layer:  p0  =  0.27 


g/cm3,  L0  =  120  pm.  The  time  interval  between  the  shaping  pulse  and  the  drive  laser 
pulse  should  be  found  from  the  requirement  that  the  two  shock  waves  merge  when  they 
hit  the  payload. 

Our  numerical  simulations  were  made  with  the  fast2d  hydrocode  developed  at 
NRL32  (more  details  about  the  code  and  further  references  are  given  in  Refs.  1,  3).  We 
ran  the  code  in  a  1-D  mode  in  order  to  verify  the  above  theoretical  considerations  and  to 
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do  the  matching.  Figure  1  shows  the  density  profiles  resulting  from  the  irradiation  of  a 
120  |tm  thick  CH  foam  layer  (p0  =  0.39  g/cnr  )  with  a  t  =  0.325  ns  long  shaping  pulse 

at  laser  wavelength  X,  =  0.248  pm  and  intensity  I p  =4.1  TW/cnr,  delivering  the 
energy  E p  =  1.33  kJ/cnr.  Propagation  of  the  shock  wave  shown  in  Fig.  1  is  well 
approximated  by  the  expression  (1)  with  a  =  0.685  and  q0  =  0.604  ,  which  can  be  re¬ 
written  as 

/  2  \l/3 

E  '  K.l/cm  n  m  si 

*,(0  =  -13x  ' - 3  x (t, ns)~° 0183 x(t,ns)°  685  pm.  (9) 

v  P o.g/cnr 

Having  varied  E  and  p0,  and  we  have  checked  that  Eq.  (9)  is  indeed  a  good  wide- 

range  approximation  for  a  CH  foam  impulsively  loaded  with  a  short  laser  pulse.  Note  that 
the  self-similarity  exponent  a  is  outside  of  the  range  in  Eq.  (3),  although  fairly  close  to  its 
upper  boundary,  which  indicates  that  the  shock  dynamics  is  mainly  determined  by  the 
laser  energy  deposited  at  the  surface,  and  is  almost  independent  of  the  pulse  duration.  The 
small  deviation  from  (3)  is  not  surprising,  since  the  simulation  deals  with  approximate 
rather  than  exact  self-similarity  of  the  flow:  the  equation  of  state  is  realistic,  not  an  ideal 
gas  model,  and  the  shock  compression  gradually  decreases  as  the  shock  wave  slows 
down;  see  Fig.  1 .  Of  course,  a  planar  blast  wave  launched  by  deposition  of  the  same 
energy  per  unit  area  would  propagate  faster,  at  the  times  of  interest,  primarily  due  to  a 

larger  value  of  the  coefficient  £o.  We  have  checked  this  by  making  the  plasma  bounded 

by  a  rigid  wall  at  x  =  0,  and  found  that  the  shock  trajectory  is  well  approximated  by  the 
formula  like  (9)  with  no  dependence  on  t,  the  appropriate  power  a  =  2/3  and  the 
coefficient  about  30  instead  of  13,  in  agreement  with  the  blast  wave  solution  (2)  of  Ref. 
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3 1 .  Deviation  from  the  exact  self-similarity,  particularly  the  gradual  reduction  of  shock 


strength  and  density  compression  ratio,  is  the  reason  why  we  chose  a  somewhat  higher 
foam  density  than  that  given  by  (8),  in  order  to  ensure  the  shock  impedance  matching 
after  the  first  shock  wave  travels  120  pm. 

Figure  2  shows  the  shock  dynamics  for  the  case  when  the  above  shaping  pulse  is 
followed  by  the  drive  pulse  after  At  =  12  ns  (the  interval  between  the  shaping  pulse  and 
the  drive  pulse  reaching  half-maximum).  For  the  drive  laser  pulse,  we  have  taken  a  4-ns 
long  Nike  pulse27  without  a  foot  [Fig.  3(a),  dashed  line].  The  density  profiles  [Fig.  2(a)] 
and  the  x-t  diagram  [Fig.  2(b)]  demonstrate  a  good  matching,  both  in  time  (the  two  shock 
waves  merge  when  they  hit  the  interface  between  the  foam  layer  and  the  payload)  and  in 
shock  impedance  (no  shock  reflection  is  observed).  In  the  next  section  we  will 
demonstrate  that  the  density  profile  has  the  desired  properties  needed  to  reduce  the  laser 
imprint. 

Fig.  3(a)  shows  the  shape  of  the  Nike  pulse  and  Fig.  3(b)  -  the  corresponding  x-t 
diagrams  for  a  62-pm  solid  plastic  target.  For  this  case  we  used  the  Nike  pulse  shape 
with  a  foot  to  control  the  target  adiabat.  We  see  that  the  ablation  front  in  Fig.  3(b)  is 
never  concave,  which  corresponds  to  its  acceleration  (in  the  negative  x-direction).  In 
contrast  with  this,  the  ablation  front  shown  in  Fig.  2(b)  is  visibly  concave  before  any 
shock  reaches  the  payload,  thus  indicating  effective  deceleration,  which  makes  the  areal 
mass  perturbations  oscillate  rather  than  grow.  Thus,  the  imprint  is  inhibited  before  the 

-5 

target  acceleration  starts,  in  a  similar  manner  as  in  graded  density  targets.  Simulations 
reported  in  the  next  Section  confirm  this  statement. 
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III.  STABILIZING  EFFECT  OF  A  DYNAMICALLY  TAILORED  DENSITY  PROFILE 


This  Section  investigates  the  efficiency  of  the  density  profile  shaping  described 
above  for  laser  imprint  mitigation.  We  now  use  the  the  fast2d  hydrocode  in  a  2-D  mode. 
The  radiation  package  of  fast2d  has  not  been  invoked  in  the  simulations  reported  here. 

In  our  previous  modeling  ’  we  have  verified  that  radiation  has  only  a  small  effect  on  the 
acceleration  of  planar  plastic  targets.  To  quantify  the  stabilizing  effect  of  density  tailoring 
with  a  shaping  laser  pulse,  one  needs  some  reference  points  for  comparison.  As  first 
shown  theoretically  and  then  confirmed  experimentally,  the  structure  of  the  sandwich 
plastic-foam  target  (see  Ref.  3,  target  2,  and  Fig.  4  therein)  is  already  much  more  imprint- 
resistant  than  a  uniform  solid  target.  All  the  targets  are  accelerated  with  a  4-ns  Nike 
pulse27  preceded  by  a  ~3-ns  low-intensity  foot:  Al  =  0.248  pm,  with  peak  intensity  50 
TW/crn2,  time  t  =  0  corresponds  to  the  instant  when  the  laser  intensity  reaches  half 
maximum.  Duration  of  the  laser  shaping  pulse  is  varied  from  0.325  ns  to  1.3  ns,  while  the 
energy  delivered  by  the  shaping  pulse  to  the  target  is  kept  constant,  1.33  kJ/cnr. 

Figure  4(a)  compares  the  perturbation  growth  in  a  sandwich  target  with  and 
without  the  shaping  pulse,  and  in  a  solid  plastic  target.  Figure  4(b)  compares  the  results 
for  a  0.325-ns,  a  1.3-ns,  and  a  0.1  ns  shaping  pulse.  The  areal  mass  variation  amplitude 
in  Fig.  4(a)  is  normalized  to  the  solid  plastic  density,  1.07  g/cm  ,  and  thus  expressed  in 
microns  (not  to  be  confused  with  the  amplitude  of  interfacial  ripples,  see  discussion  in 
Ref.  24  for  details).  For  the  simulations  in  Fig.  4  the  laser  nonuniformity  is  modeled  as  a 
0.2%  single-mode,  constant-phase  sine-wave  intensity  perturbation  at  A  =  30  pm 
superimposed  on  an  otherwise  spatially  uniform  drive  laser  pulse.  The  intensity 
perturbation  in  the  shaping  pulse  is  introduced  in  a  similar  way,  and  the  curves  in  Fig. 
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4(a)  show  the  effect  of  varying  the  shaping  pulse  perturbation  amplitude  from  0  to  0.2% 
to  2%.  Figure  4(a)  also  presents  the  corresponding  time  histories  obtained  with  the  same 
drive  pulse  (but  without  any  shaping  pulse),  for  the  same  plastic-foam  sandwich  target 
and  for  a  62 -pm  solid  plastic  target. 

The  solid  plastic  target  is  strongly  distorted  by  the  drive  pulse.  If  the  drive  pulse 
lasted  any  longer,  a  bubble  would  have  punctured  through  the  target.  The  foam-plastic 
sandwich  target  without  the  shaping  pulse  is  somewhat  less  distorted.  However,  the  areal 
mass  oscillations  produced  by  the  drive  pulse3  are  so  strong  that  this  target  —  with  a 
relative  mass  variation  slightly  below  20%  -  is  not  likely  to  survive  much  more 
acceleration.  A  shaping  pulse  with  a  2%  intensity  variation  brings  a  considerable 
improvement  —  by  a  factor  of  2-3.  Not  unexpectedljthe  best  results  are  found  for  the 
0.2%  and  0%  intensity  variation  in  the  shaping  pulse.  Surprisingly,  the  mass  variation 
amplitude  in  the  case  of  a  0.2%  is  even  less  than  it  would  be  with  a  perfectly  smooth 
shaping  pulse:  0%!  The  explanation  is  simple:  with  the  intensity  variation  given  by  the 
same  sine  wave  function  for  both  the  shaping  pulse  and  the  drive  pulse,  the 
hydrodynamic  perturbations  caused  by  these  pulses  are  coherent,  they  can  interfere,  and 
in  this  case  they  are  seen  to  interfere  negatively. 

Shorter  shaping  pulses  in  Fig.  4(b)  provide  better  stabilization.  This  is  explained 
by  the  properties  of  the  mass  perturbation  growth  caused  by  a  single-mode  laser  intensity 
modulation  during  the  shock  transit  time,  which  have  been  established  in  the  studies  of 
laser  imprint.  Let  us  assume  that  the  shaping  pulse  duration  r  is  shorter  than  the 
decoupling  time  tic  between  the  critical  surface  and  the  ablation  front:  tdc  >  r .  This 
decoupling  time  is  the  characteristic  smoothing  time  of  the  pressure  modulation  at  the 
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ablation  front  caused  by  the  laser  intensity  modulation  in  the  absorption  zone.  The 
smoothing  of  pressure  modulation  is  due  to  the  buildup  of  the  plasma  corona  separating 
the  two  regions,  i.  e.  tdc  =  \/kvc,  where  k  is  the  perturbation  wave  number,  and  vc  is  the 
velocity  of  the  absorption  zone  with  respect  to  the  ablation  front.  Then  the  time 
dependence  of  the  areal  mass  modulation  amplitude  could  be  presented  as22 
o  SI 

Sm(t)  =  — x — —xF(kUt) ,  (TO) 

k  :p 

where  81  II  is  the  relative  modulation  amplitude  of  the  shaping  beam  intensity,  U  is 
the  mass  velocity  of  the  shocked  plasma,  F(  t )  is  a  function  of  normalized  time  t  =  kUt  . 
For  the  sake  of  simplicity  we  neglect  the  weak  parametric  dependence  of  F(t )  on  the 
laser  intensity  via  the  blowoff-plasma-to-ablation- front  density  ratio.  Then  for  the  short 
times  of  interest,  F(  t )  is  simply  a  growing  function  of  its  dimensionless  argument, 
approximately  proportional  to  7 2  ■  The  mass  velocity  U  scales  as  p'a2 ,  where  pa  is  the 
ablative  pressure  driving  the  shock  wave,  which,  in  turn,  scales  as  some  power  of  the 
intensity  /  of  the  shaping  laser  beam:  Pa  x  I], ,  with  v  somewhere  between  0.66  and 

0.8.  We  compare  the  shaping  pulses  at  constant  energy  E  per  unit  area, 

=  Ep  It  °<t  1 ,  so  that  jj  cc  x  vil .  The  end  of  the  shaping  pulse  t  =  x  thus  corresponds 
to  a  dimensionless  time  tp=kU'C  ' /2 .  Since  V  /  2  <  1 ,  Eq.  ( 1 0)  demonstrates  that  a 
shorter  shaping  pulse  will  cause  smaller  perturbation  imprint  into  the  flow  at  the  end  of 
the  pulse  Sm(t  =  z)oc  r2  v .  Moreover,  even  the  time  derivative  of  the  mass  perturbation  at 
the  end  of  the  shaping  pulse  is  smaller  for  a  shorter  pulse.  This  conclusion  does  not 
change  if  we  take  into  account  the  weak  parametric  dependence  of  F(  t )  on  the  laser 
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intensity  via  the  density  ratio  and  the  decoupling  of  the  intensity  modulation  in  the 
shaping  beam  from  the  pressure  perturbations  at  the  ablation  front.  Both  effects  are  more 
pronounced  at  higher  laser  intensities  and  tend  to  decrease  the  mass  variation. 

Realistically,  a  target  cannot  be  manufactured  perfectly  smooth.  Since  the  surface 
finish  and  the  intensity  variation  in  the  beam  are  both  relatively  small  and  uncorrelated, 
their  effects  could  be  studied  separately.  Figure  5  demonstrates  the  perturbation  growth 
caused  by  placing  a  single-mode  ripple,  A  =  30pm,  peak-to-valley  amplitude  2 a0  =  2  pm 
on  the  foam  surface  (an  areal  mass  modulation  equivalent  to  0.79  pm  in  solid  CH),  while 
the  shaping  pulse  and  the  drive  beam  are  held  perfectly  smooth.  A  2  m  perturbation  is 
much  worse  than  the  surface  finish  expected  for  an  ICF  target.  We  observe  that  the 
perturbation  with  the  laser  shaping  pulse  grows  to  almost  as  large  an  amplitude  as 
without  it.  In  this  case  the  shaping  pulse  is  seen  to  cause  strong  sonic  oscillations  in  the 
shocked  plasma,  which  amplify  the  mass  variation  in  the  ripples  by  a  factor  of  3-4  and 
persist  until  the  drive  laser  pulse  arrives  to  provide  a  higher  (in  comparison  with  the 
original  surface  ripples)  initial  perturbation  for  the  subsequent  growth.  This  rather 
unexpected  behavior  is  because  the  perturbed  flow  contains  a  rippled  expansion  wave, 
which  immediately  follows  the  rippled  shock  wave.  Oscillations  of  a  rippled  shock  wave 
are  known  to  decay  from  their  peak  initial  amplitude  to  the  lower  values,  see24  and 
references  therein.  On  the  contrary,  some  amplification  —  by  a  factor  of  three  or  so  -  of 
the  initial  mass  variation  amplitude  is  characteristic  of  a  rippled  rarefaction  wave.  ’  The 
large-amplitude  oscillations  seen  in  Fig.  5  after  the  shaping  pulse  are  apparently  caused 
by  the  expansion  wave  component  of  the  flow.  From  our  simulations  we  conclude  that 
the  rms  surface  finish  at  A  =  30  |im  should  not  exceed  0. 1  pm,  which  is  consistent  with 
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the  conventional  smoothness  requirements  for  direct  drive  laser  fusion  targets.  Note  that 
no  additional  requirements  are  made  for  the  rear  (inner)  surface  roughness,  as  the  shock 
wave  driven  by  the  drive  laser  pulse  merges  with  that  driven  by  the  shaping  pulse  before 
it  hits  the  rear  surface. 

Now  consider  the  residual  non-uniformity  of  both  the  shaping  pulse  and  the  drive 
laser  pulse  left  with  ISI  smoothing.  The  theoretical  model  of  this  smoothing  technique 
and  its  numerical  implementation  in  fast2d  are  described  in  detail  in  Ref.  26.  We  assume 
a  bandwidth  of  1  TW,  which  implies  a  mis  intensity  variation  in  a  single  Nike  beam  time- 
averaged  over  4  ns  of  less  than  2%.  With  N  B  statistically  independent  overlapping 
beams,  the  intensity  variation  in  wavelengths  greater  than  10  m  decreases  further  by  a 
factor  of  Nb112  ,  to  about  0.25%  for  Nike  s  NB  =  36  .  For  the  shaping  pulse,  which  is 
shorter  than  4  ns,  the  time-averaged  rms  intensity  variation  in  a  single  beam  is  higher  by 
a  factor  of  (4  ns/r)1  2  ~  2  to  3.5. 

As  far  as  the  energy  is  concerned,  one  beam  out  of  36  (with  a  pulse  shortened  to 
1  ns  or  less)  would  be  sufficient  to  provide  the  desired  energy  of  the  shaping  pulse,  since 
the  shaping  pulse  needs  less  than  1%  of  the  laser  energy  in  the  drive  pulse.  However,  a 
single  beam  is  6  times  less  uniform  than  36  overlapping  beams.  The  question  is  whether  a 
single  beam  is  sufficiently  smooth  not  to  pre-introduce  the  imprint  before  the  drive  pulse. 

Figure  6  answers  this  question.  It  compares  the  cases  where  [Fig.  6(a)]  the 
shaping  pulse  is  delivered  by  a  single  beam,  ISI-smoothed  at  1  THz  to  cases  [Fig.  6(b)] 
with  36  such  beams,  overlapped  and  carrying  the  same  energy  to  the  target.  The  drive 
laser  beam  in  both  cases  consists  of  36  beams.  We  see  that  a  single  beam  used  as  a 
shaping  pulse  is  definitely  not  smooth  enough  to  mitigate  the  imprint.  The  single  beam 
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imprints  enough  perturbations  of  its  own,  that  the  resulting  mass  variation  amplitude  is 
even  larger  than  in  the  absence  of  a  shaping  pulse.  On  the  other  hand,  with  the  shaping 
pulse  sufficiently  smoothed  (  that  is,  with  the  same  number  of  beams  as  in  the  drive 
pulse),  it  obviously  helps  to  mitigate  the  imprint  caused  by  the  non-uniformity  in  the 
drive  pulse,  Fig.  6(b).  Reduction  of  the  non-uniformity  of  a  single  laser  beam  by  almost 
another  order  of  magnitude  through  increased  bandwidth  is  hardly  practical,  see  Ref.  24. 
Therefore  the  shaping  pulse  should  be  formed  by  overlapping  multiple  optically 
smoothed  laser  beams. 

III.  CONCLUSIONS 

We  have  demonstrated  that  the  mass  imprint  can  be  reduced  by  use  of  a  single 
shaping  laser  pulse,  with  the  effect  as  large  as  is  predicted  for  a  tailored  density 
profile.  The  energy  of  the  shaping  pulse  can  be  a  small  fraction,  less  than  1%,  of  the  total 
energy  in  the  drive  laser  beam.  The  shaping  pulse  must  be  short  on  the  nanosecond 
time  scale:  1  ns  duration  will  do,  0.3  ns  is  better. 

We  found  that  the  shaping  pulse  must  also  be  smoothed  sufficiently  to  mitigate 
the  imprint.  Otherwise,  the  shaping  pulse  imprints  its  own  non-uniformity,  which  persists 
during  the  time  interval  between  the  shaping  pulse  and  the  drive  pulse  in  the  form  of 
strong  sonic  waves,  providing  seeds  for  the  perturbation  growth  associated  with  the  drive 
pulse.  Similarly,  the  shaping  pulse  can  increase  the  initial  surface  roughness  in  the  target, 
transforming  this  initial  mass  variation  into  strong  sonic  waves  on  the  target  surface. 
There  is  a  corresponding  upper  limit  on  the  surface  finish  of  the  target.  The  requirements 
on  the  smoothness  of  the  shaping  pulse  and  the  target  surface  are  strict,  but  not 
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unrealistic.  In  particular,  the  laser  drive  uniformity  available  with  36  overlapping  Nike 
beams,  ISI-smoothed  at  1  THz,  seems  to  be  sufficient  for  this  purpose. 

High-gain  or  ignition  direct-drive  targets  that  use  this  method  for  mitigating  the 
laser  imprint  could  be  designed,  but  this  task  is  beyond  the  scope  of  the  present  paper.  A 
necessary  intermediate  step  is  designing  and  performing  an  experiment  to  demonstrate  in 
planar  geometry  the  efficiency  of  this  method  of  imprint  mitigation.  This  work  is  now  in 
progress.  It  is  too  early  to  extrapolate  to  large-scale  facilities,  and  we  limit  ourselves  to 
the  following  observation.  Since  the  energy  in  the  shaping  pulse  is  less  than  1%  of  the 
energy  in  the  drive  beam,  our  method  of  stabilization  at  the  ignition  energy  level  (~  1  MJ 
in  the  drive  pulse)  would  require  something  like  a  36-beam,  10-kJ  laser,  somewhat  more 
powerful  than  Nike  but  a  lot  less  powerful,  less  expensive  and  less  complicated  than  the 
~100  kJ  picosecond  laser  needed  to  pursue  the  fast  ignition  option.34 
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Figure  captions 


Fig.  1.  Time  evolution  of  the  density  profiles  in  a  foam-plastic  sandwich  target  irradiated 
by  a  0.325  ns  long  shaping  pulse  at  t  =  -12  ns.  The  profiles  are  shown  with  a  1  ns  interval. 
Fig.  2.  Time  evolution  of  the  density  profiles  (a)  and  the  x-t  diagram  showing  shock 
dynamics  (b)  for  the  same  conditions  as  in  Fig.  1  but  with  the  drive  laser  pulse  following 
the  shaping  pulse.  The  upper  envelope  curve  of  the  high-density  area  is  the  ablation  front. 
When  the  shaping  pulse  is  applied,  this  curve  is  concave  at  early  time,  which  implies 
deceleration. 

Fig.  3.  (a)  Shape  of  the  standard  Nike  pulse  with  a  low-intensity  foot  (solid  line)  and 
without  a  foot  (dashed  line,  shifted  by  ~0.5  ns  closer  to  the  shaping  pulse,  to  provide 
matching  under  the  conditions  of  Fig.  1 ).  /  =  0  is  the  moment  when  the  Nike  pulse  with  a 
foot  reaches  half-maximum.  The  shaping  pulse  starts  at  t  =  -12  ns.  (b)  The  x-t  diagram 
showing  shock  dynamics  for  a  62  pin  solid  plastic  target  irradiated  with  a  Nike  pulse  and 
no  shaping  pulse.  In  this  case  the  ablation  front  is  never  concave. 

Fig.  4.  (a)  Evolution  of  mass  variation  amplitude  in  a  sandwich  target  with  and  without 
the  0.325  ns  shaping  pulse  (the  single-mode,  constant  phase  intensity  variation  is  the 
same  in  the  shaping  pulse  and  the  drive  pulse)  and  in  a  solid  plastic  target  without  a 
shaping  pulse,  (b)  Same  comparing  the  effects  of  a  0.325  ns,  1.3  ns  and  0.1  ns  shaping 
pulses  with  a  2%  intensity  modulation  in  a  sandwich  target.  All  the  other  simulations  in 
this  paper  utilized  the  0.325  ns  shaping  pulse. 

Fig.  5.  Evolution  of  mass  variation  amplitude  in  a  sandwich  target  rippled  from  the 
surface  with  and  without  the  0.325  ns  shaping  pulse;  both  the  shaping  pulse  and  the  drive 
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pulse  are  perfectly  uniform.  The  shaping  pulse  is  seen  to  produce  strong  sonic  oscillations 
whose  amplitude  exceeds  the  initial  mass  variation  due  to  the  surface  ripples. 

Fig.  6.  Evolution  of  mass  variation  amplitude  in  a  sandwich  target  with  both  shaping  and 
drive  laser  pulses  ISI  smoothed  at  1  THz  bandwidth,  (a)  The  shaping  pulse  is  a  single 
beam,  the  drive  pulse  is  36  overlapping  beams,  (b)  Both  the  shaping  and  the  drive  pulses 
are  36  overlapping  beams,  (c)  No  shaping  pulse. 
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Metzler  et  al..  Fig.  2(a) 
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Metzler  et  al..  Fig.  4(a) 
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Fig.  4(b). 


Metzler  et  al..  Fig.  4(b) 
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Metzler  et  al..  Fig.  5 


